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Abstract 

We compute the final kinetic energies of the fragments emitted in the 
light charged particle accompanied cold fission of 252 Cf taking into ac¬ 
count the deformation and the finite-size effects of the fragments and 
integrating the equations of motion for a three-body system subjected 
only to Coulomb forces. The initial conditions for the trajectory cal¬ 
culations were derived in the frame of a deformed cluster model which 
includes also the effect due to the absorbative nuclear part. Although 
the distributions of initial kinetic energies is rather broad we show that 
in cold fission the initial conditions can be better determined than in 

the usual spontaneous fission. 
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I. INTRODUCTION 


The cold ternary fission is a rearrangement process of a large group of nucleons 
from the ground state of the initial nucleus to the ground state of the three final 
fragments. Like in the case of spontaneous and thermal-induced fission a ternary 
component of a few tenths of percent is present also in the cold fission process [|T|-j3| . 

In order to determine the configuration and the dynamics of the fissioning nu¬ 
cleus at scission, the experimental data for the light charged particle (LCP) emitted 
in the fragmentation process are analyzed and compared with the theoretical results 
obtained via trajectory calculations. In the past a large number of studies were 
devoted to the trajectory calculation, specially for a-particles in the point charge 
approximation and without the account of nuclear forces. The alphas were consid¬ 
ered to be emitted from the neck region [ffj-pj. There have been also some authors 
who considered the finite size and the deformation effects EMI and showed that 
these geometrical factors are influencing sensitively the angular distributions of the 
LCP. 

However in all these approaches to the ternary spontaneous fission the problem 
of choosing the initial parameters of the trajectory calculations is complicated by the 
fact that various theories give different predictions. Since only the energies and the 
angles of the three particles can be experimentally determined, solving the equations 
of motion backward in time will not provide a full information on the geometric and 
dynamic characteristics of the fissioning system at the moment of the LCP-emission. 
The only possibility is to probe various combinations of assumed initial conditions 
and then compute the trajectories for comparison with the available experimental 
data. In the hot ternary fission the initial conditions are so numerous that in or¬ 
der to encompass as much as possible combinations, Monte-Carlo techniques were 
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employed ||TT| . p~2 


For the cold ternary fission the initial conditions are better known |3].|T4|. First of 
all the fragment deformations are those of their ground states. This fact prompted 
us to calculate the final characteristics of the LCP emitted in the cold fission of 
252 Cf for different mass splittings, and see how the static deformations and the 
finite size are modifying the outcome of the trajectory calculation. Using forces, 
computed through a double folded integration of the Coulomb interaction between 
two quadrupole deformed heavy ions, we derived the equations of motion for the 
three-body problem, and solved them numerically . The solution of this set of 
equations provided the final angle of the LCP with respect to the fission axis and 
its kinetic energy. We compared our calculations with the point-like trajectory and 
some experimental consequences were discussed. 


II. DETERMINATION OF THE INITIAL CONDITIONS 

Usually, in trajectory calculations for the spontaneous fission different choices 
are taken for the initial kinetic energies of the fragments emitted in the process. 
For example the initial kinetic energy of the two main fragments and of the a 
emitted in the spontaneous ternary fission should be around 0.5 MeV according to 
the statistical theory and the equipartition principle [|7]||. On the contrary, in the 
dynamical theory of fission |U| the nascent fragments at scission are predicted to 
be moving with appreciable kinetic energy (20-50 MeV). The initial velocities of 
the heavy fragments are considered to have non-zero components only along the x- 
axis. The initial velocity of the light fragment Ul( 0) is related to the initial velocity 
of the heavy fragment u#(0) in such a way that the total momentum of the two 
fission fragments is zero along the x-axis, i.e. Vl(0) = j^Vh( 0). This reasonable 
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assumption will be applied also by us. 

Therefore we have to determine the following initial conditions : a) The tip 
distance d\ b) the kinetic energies of the two main fragments Eho, El 0 and the 
kinetic energy of the LCP Elcp ; c) the initial geometric configuration of the LCP, 
i.e. the position between the fragments and the angle d° LCP between the direction 
of motion of the LCP and the axis joining the two main fragments; d) The shape of 
the fragments. 

The determination of these quantities in the ternary cold fission will be facili¬ 
tated, up to a certain extent, by the peculiar characteristic of the process, i.e. the 
fragments are emitted with total kinetic energy TKE close to the corresponding 
ternary decay energy Q t . In order to achieve such large TKE values, the three final 
fragments should have very compact shapes at the scission point and deformations 
close to those of their ground states, similar to the case of the cold binary frag¬ 
mentations. One may next suppose that the shapes of the fragments will not be 
modified when the fragments move away in the Coulomb field of each other. Thus, 
the problem of the fragments shape in the initial configuration is easily established 
for the cold fission. 

In order to determine the kinetic energies of the two main fragments we make 
use of the considerations derived from the deformed cluster model that we employed 
in previous papers for the study of the ternary cold fission 


. In this deformed 


cluster model the barrier between heavy fragments (for binary fission) and the barrier 
between the LCP and the heavier fragments (for ternary fission) can be calculated 
quite accurately due to the fact that the touching configurations are situated inside 
of the barriers. For the two fragments, the exit point from the potential barrier is at a 
tip distance d around 3 fm, as can be seen in Figure 1, for the case 248 Cm —> 104 Mo 
+ 144 Xe. This barrier is much thiner than the barrier between the LCP and the 
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heavier fragments, and thus in our model first the two heavier fragments penetrate 
the potential barrier between them and later on the LCP is emitted. Consequently 
the mass distributions of the heavier fragments are very similar to those of the cold 
binary fission of an initial nucleus leading to the same heavy fragments, i.e. 248 Cm 


if the LCP is an a 0- or 242 Pu if the LCP is 10 Be [|nj. The decay energy for such 
a binary fragmentation will be Qlh = Qt ~ Qlcp, where Q t is the ternary decay 
energy of 252 Cf and Qlcp is 6.22 MeV for a and 8.71 MeV for 10 Be. 

On ground of the cold fission characteristics mentioned above one may conjecture 
that at the exit point (second turning point) of the two heavier fragments, their 
potential energy is equal to Qlh and their kinetic energy is equal to zero. When the 
fragments move apart, i.e. their tip distance increase, their kinetic energy increase 
to. In order to estimate the total kinetic energy of the fragments we have to find 
out at which tip distance the release of the LCP is likely to occur and compute at 
that point the potential energy, i.e. 


TKE(d) = TKE l + TKE h = Q LH - V LH (d) 


( 1 ) 


Using the conservation of linear momentum invoked above we have 

A h 

TKE l = -t—TKE h (2) 

Al 

and the individual kinetic energies in terms of the total kinetic energy read 

TKE r = t 1' , TKE(d) (i = L, H) (3) 

Now we turn to the problem of determining the tip distance d. It is reasonable to 
suppose that d should correspond to the configuration at which the LCP is released. 
In Figures 2a and 2b we plotted the ternary potential seen by the LCP (in this case 
an a) in the field of the two heavy fragments. As we shall se later the LCP should 
stay between the two heavy fragments in a position which should avoid its absorbtion 
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by any of the fragments. We see in Figure 2a, that for tip distances up to 7 fin, the 
a is facing a thick barrier in the transversal direction. Eventually as the distance 
between the fragments increases the pocket in which the a is located becomes more 
and more shallower untill it disapears around d = 8 fm. Therefore one may conclude 
from these qualitative arguments that the initial tip distance between the two main 
fragments should not be larger than that corresponding to the disapearence of the 
LCP pocket. On the other hand for tip distances smaller than 6 fm the emission 
of the a is strongly hindered by a thick barrier even for a rather high zero energy 
E° > 3MeV (see Fig.3). 

If we choose d = 8 fm for the example considered in Figure 1, then we get for 
the total kinetic energy of the two main fragments TA'A=46.21 MeV which is much 
larger than the corresponding kinetic energy in the spontaneous fission. For d — 6 
fm the kinetic energy will drop to TKE— 28.78 MeV. Repeating this calculation for 
other mass splittings we conclude that the kinetic energy of the main fragments is 
ranging in the broad interval 25 - 50 MeV, but as we shall see bellow it is correlated 
to the kinetic energy of the emitted alfa particle through the tip distance. 

We are left now with the determination of the LCP geometrical and dynamical 
initial characteristics. For that we invoke a receipt proposed by Boneh et al. || which 
consider as a possible choice for the LCP position, the point of minimum potential 
energy (the saddle point of the potential energy surface). If the heavy-fragments 
would have to interact with the LCP via point-like Coulomb forces, this electrostatic 
saddle point would be determined by Z H /R 2 aH = Z L /R 2 aL) where R a i{i = L,H) is 
the distance between the LCP and the main fragment i. It is readily seen from 
Figure 2a,b that in the case of our deformation dependent cluster model, where the 
nuclear forces are introduced via the M3Y potential, this saddle point corresponds to 
the position where the combined Coulomb and nuclear forces exerted by the heavy 


6 


fragments on the LCP cancel each other and the potential surface will have a relative 
minima at this point. To establish more precisely the location of this electro-nuclear 
saddle point, we use the multipolar decomposition for the M3Y potential [|T6[, and 
the above condition translates to 


v dv$ = v dv™(R aL ) 
^ dR n - ^ 


\ ^ dR a L 

which is a generalization of the point-like Coulomb equilibrum condition. In the 
laboratory frame of reference, we choose the z-axis as the initial fissioning axis of 
the two heavier fragments, with the origin at the tip of the left (heavy) fragment. 
Then the location of the electrostatic saddle point is given analitically by the formula 


( 4 ) 


/ n d + Ol + CtJl . . 

Z a {d) = - 7 - a L (5) 

1 + <Jz L /z H 

where cq (i = L , H) are the major axes of the quadrupole deformed main fragments. 
For the pair considered in Fig.1-3, z a (d) ^0.58 d using the point-like Coulomb forces 
and 0.51d in our model where nuclear forces are included to. 

As we already noted above the potential energy of the LCP positioned at the 
electro-nuclear saddle point will have a minimum in the ^-direction. It is clear that 
the LCP can have no component of its velocity along the x-axis since this would 
result in a possible absorbtion of the LCP by the deep potential wells of the two 
heavier fragments instead of being emitted []. The only possibility for the LCP to 
survive the descent of the decaying system from scission to the release point is to 
have a momenta directed only along the y-axis. As can be infered from Figure 
2 the locus of the saddle point is on the bottom of the potential well. Taking 


1 In the case we would employ forces with repulsive nuclear core the LCP will be once 
again prevented to move in the ^-direction. 
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sections of the potential surface along the y -axis at z corresponding to the saddle 
point, the resulting potential slice will look similar to a one-dimensional harmonic 
potential well (see Figure 3). When the tip distance increases, the well becomes more 


and more shallower untill it vanishes completely. Following an idea from [11 


we 


will approximate the potential Vlcp with an harmonic potential in the ^/-direction, 
centered at the saddle-point 


Vlcp ~ V sa ddie + ~jCy 


( 6 ) 


where V 3a m e - Vlcp(z - ^saddle, V = 0) and C - - n is the st.ifness. It can be 

u y y=o 

shown after some algebra that the elastic constant value is given by the expression 

A(A + 1) V\q\ (Rai)\ 


£ _ \ 1 \ ( 9V\q\ (Rai) 


R°- 

i=L,H A>0 


DR ni 


Rrvi 


( 7 ) 


where R a i is the distance from the fragment i to the LCP (a) on the z-axis : 


R 


o _ 

OlL 



,R 


o _ 

aH ~ 



( 8 ) 


where D is the inter-fragment distance. From here we get an estimation for the initial 
kinetic energy of the LCP supposing that it can be identified with the zero-energy 
in the harmonic potential well, i.e. 

Elcp = --— (9) 

2 V m LCP 

Consequently, a degree of uncertainty in the initial kinetic energy occurs also for 
the LCP. For increasing tip distance the kinetic energy of the LCP decreases. One 
might suppose that in the range 6-8 fm, for the tip distance, the LCP has the 
possibility to escape by tunneling or by the dissapearence of the barrier. Further 
the velocity corresponding to this kinetic energy, v a = v/ rra^cp l iave a nonzero 
component only with respect to the y- axis, according to the above discussion. 
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III. TRAJECTORY EQUATIONS 


In order to write down the equations of motion we have to establish the geometry 
of the system not only at the beginning but also long after the release of the LCP. The 
forces being central, and the initial velocities are in-plane the problem is simplified 
by a two-dimensional approximation. There will be required six coordinates and six 
velocities, which are governed by a system of twelve first order ordinary differential 
equations. In Fig. 4, we see the three fragments and the forces acting between them, 
just after the release of the LCP. Contrary to other works we take into account the 
forces exerted by the LCP on the fragments. After deriving the initial conditions in 
the previous section, taking into account the nuclear forces in the calculation of the 
barriers, we proceed now to the calculation of the trajectories by considering only 
the Coulomb forces. Since the kinetic energies of the fragments are rather high, this 
approximation is good even in the point-charge approximation. 


In previous papers ||,|l4],[r7|] we used a double folding potential for the heavy-ion 
interaction. Presently we shall consider only the Coulomb part of this interaction, 
between two ions, i.e. 



Pi(n)p2 

rq + R- r 2 



( 10 ) 


where P\(2){r) are the charge ground-state one-body densities of the fragments. The 
one-body densities are taken as Fermi distributions in the intrinsic frame for axial- 
symmetric nuclei 


p(r) = —(ii) 

1 + e « 

with R(9) = i?o(l + J2x>2 P\Y\o(6, 0)). In what follows we consider that the sym¬ 
metry axes of the fragments are lying in the same plane. Using the formalism 
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presented in |1(|, the interaction between two heavy ions with orientation U\,U 2 of 
their intrinsic symmetry axes with respect to the fixed frame, reads : 


V(Ri2] = £ V^(fl 12 )it I (co8Wi)i^' , (co 8 « J )P J ,(co S n 12 ) (12) 

Al, A 2 , A 3 

where 



with (coscji) and P^ M (coso; 2 ) being the associated Legendre polynomials which 
describe the relative orientation of the two fragments whereas PA 3 (cosf2 12 ) describes 
the orientation of the axis joining the two nuclei with respect to the laboratory 
frame. In the present study the LCP is spherical and thus the interaction between 
the LCP and one heavy fragment i(=L , H) will get a simplified form 

V(R ai ) = J2 V AOA {Rai)P\( cos 6 i) (14) 

A 

The following approximation can be applied for the two heavy fragments : Since 
their relative orientation does not change significantly at the beginning of the qua- 
siclassical motion, one can neglect the relative orientation of the heavy fragments: 


V(Rlh)= £ vIIX{R lh)Px z (*> S <) l „) (15) 

A 1 A 2 A 3 

The force acting between a pair of fragments can be written: 


F,j = -W(fly) 


(16) 


The force acting between the two heavy fragments is given by: 

'dV°£%( R ^) p 2 ,, , E UX(Rlh) 

-P\ 3 (cos<p) cos (p — 


Flh= —e x T, 

A 1 A 2 A 3 


dlli,n 


Go, 


^ E 

^ 1 A 2 A 3 


' dVSXlsiRLH) 

9Rlh 


P\ 3 (cos cj)) sin 4> + 


Rlh 

V&VURlh) 

Rlh 


Pl 3 (cos 0) sin 4> 


P\ 3 (cos (/)) cos (p ) (17) 
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whereas the forces exerted by the fragments on the LCP read: 
F Ha= —e x 


Y^ \ o \(RaH) D { , N ! Y^ ^AO x(Rqh) D l / , ^ . , 

V-^T77-PAfcos-i/’iJcosV’i - > „ -P A (cospi) sinpi 

A>0 9RaH A>2 R <* H 


— e 


^ 8y > o >( fl °«) p A ( cosV , 1 ) sinV , 1 + ^ y (Jg ° g) Pj(cosVi)cosV, (18) 

A>0 9RaH X>2 RaH 


F La *2 a 


Y^ a 0 A {Fql) d [ ! ^ t y^ ^AOA (Fol) d i^ ■ , 

2^ - -Pa (cos 1p2) COS 1p2 - 2^ ^ - P\[COS'tp 2 ) Slllp2 

^A>0 A>2 


— e 


Y ^A 0 0° D ( i \ ■/ , Y^ ^A 0 \( R cil) n l/■ , \ , 

y | > -yb- P A (cos + V-^- P A (cos ip 2 ) cos tp 2 

A>0 A>2 RaL 

The equations of motion of the three nuclei are: 


(19) 


MlVl = F lh - F La 
M H r H = -F lh - F Ha 
W^aVa. FL a T Fj{ a 


( 20 ) 


Here we assumed that the two heavy fragments have the same multipolarity in 
deformations. In this paper we consider only quadrupole deformations. 

The above system was solved numerically employing the lsoda package for or¬ 
dinary differential equations, with automatic method switching for stiff and nonstiff 
problems [18| . 

In Figure 5 a, b, c we presented the trajectory of the three fragments for the two 
extreme initial conditions (with high and with low kinetic energies of the heavier 
fragments) in a sequence of 10 time steps. The time scale is divided into increments 
of At = 1.8 x 10 22 . In figure 5a we display the trajectories of one of the most 
asymmetric splittings, recorded in experiment, i.e. 152 Nd + 92 Kr. Since in this case 
the a feels a stronger repulsion from the heavy fragment, it will be deflected at a 
larger angle in the direction of the light fragment. In the case of the splittting 144 Xe 
+ 104 Nd this deflection will be less pronounced (Fig.5b) and for the more equilibrated 
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splitting, i.e. 132 Nd + 116 Pd, the a will be only slightly deflected (Figure 5c). We 
thus observe that in all the cases the a-particle is deflected in the direction of the 
light fragment, but with a larger angle when the initial kinetic energy of the heavier 
fragments is higher. This fact should be attributed to the low energy of the a (< 
IMeV) which makes it to feel for a longer time the repulsion coming from the heavy 
fragment. In Table I we present the final kinetic energy Ef and the asymptotic 
angle 6£ for the three splittings mentioned above when we employ point-like and 
size dependent forces. In all cases we observe the decreasing of Ef with increasing 
tip distance d. The mean of these two energies is not far from the value of ( Ef) 
= 16.0 ±0.2 MeV which is the most probable a-particle energy predicted by the 
trajectory calculation for the hot fission. Thus, the phenomenon of a-particle energy 
amplification in the cold fission seems to follow the same pattern like in normal 
fission. This effect should be attributed solely to the predominant effect of the 
Coulomb field and less to deformation or finite size effects. It should also be remarked 
the near constancy of the final LCP kinetic energy for different mass splitings at the 
same tip distance, a fact already remarked long time ago in spontaneous fission [BJ. 
In what concerns the angles at which a-particles are emitted their dependence on 
the mass splitting is obvious. Deviations from the axis perpendicular to the fission 
axis increase with the mass ratio. The difference observed between the two sets of 
data points to an important influence of the geometrical factors, which however does 
not alter the general trends of the process. 

Naturally one might next ask if the experimental status of the problem allows 
the comparison with the results presented in this paper. Up to now there are no 
special data on cold fission available since from the experimental side it would mean 
to set a trigger on neutronless events, which is very difficult to attain in practice. 
There are available data on the alpha (and other particles) spectra, as a function 
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of the total excitation energy, reaching TXE = 10 MeV within the experimental 
accuracy performed by the Darmstadt group with the DIOGENES setup, and in a 
more recent work at the MPI Heidelberg jl9|. These data does not contain special 
effects in the alpha spectra, when the cold fission regime is approached, except that 
the mean energy increases nearly linearly with decreasing TXE. This would mean 
that if the linear dependency would be extrapolated to TXE = 0 MeV, i.e. when 
the scission configuration tends to become compact, like in our model, the average 


kinetic energy of the a will approach the value 17.5 MeV [0. According to the 
calculations presented in this work a range between 12 MeV to 20 MeV should 
be expected for the final kinetic energy if we consider that the a particle occupies 
the lowest states in the pocket formed from the interaction with the two heavier 
fragments. Therefore the experiment doesn’t show a distinctive a kinetic energy 
distribution for cold fission, a fact which is in agreement with the calculations we 
presented in this paper. In order to establish more precisely ( E a ) we should carry 
out Monte-Carlo calculations. The fact that the experimental value is slightly higher 
than in hot fission (15.9 MeV) is a sign that the a is emitted earlier in cold fission, 
according to the uncertainty relation for energy A E ■ At ~ fi. 


IV. FINAL REMARKS 

We presented a receipt to determine the initial conditions for trajectory calcula¬ 
tions in the ternary cold fission of 252 Cf. Compared to the case when the fragments 
are emitted with high excitation energy, the initial conditions are restricted to a 
narrower range of values as a consequence of the pecularities of the process, mainly 
the compact shape of the fragments. 

In our model the a-particle cannot be emitted at a tip distance larger than 8 


13 










fm, because as we showed for such a distance the a particle is no longer under the 
influence of the attractive nuclear forces, and on the other hand we disregard tip 
distances smaller than 6 fm because the LCP wavepacket hlling the lowest state in 
the potential well has a small probability to tunnel through the thick transversal 
barrier. 

The location of the LCP was fixed at the electro-nuclear saddle point, which is 
also model dependent. Due to the finite size and the deformations of the fragments 
this location will be shifted with respect to the location of the electrostatic saddle 
point, and consequently the outcome of the trajectory calculation will be altered to 
a certain extent. 

The initial kinetic energy of the LCP was considered to coincide with the lowest 
level occupied in a one-dimensional harmonic potential well oriented perependiculary 
to the fission axis. This energy is decreasing with the tip distance. 

As have been pointed by Halpern , there is no reason to believe that the third- 
particle ejection rates should be independent of the initial angular momentum. In 
our case the spin of the parent nucleus ( 252 Cf) being zero the angular momentum 
imparted to the fragments and their relative angular momentum is mainly due to 
the creation of a molecular configuration at the scission point [^T]]. In the model 
presented in this paper we didn’t took into account the influence of collective molec¬ 
ular excitations, like bending or wriggling, nor the torques exerted between the 
fragments during the quasiclassical motion. The inclusions of these suplementary 
degrees of freedom could alter the initial configuration. This would be an interesting 
topic for a future investigation of the scission configuration in ternary cold fission. 
Moreover the evolution from scission to the release point of the LCP should be done 
in a dynamical way, i.e. writing equations of motion not only for the translational 
and rotational degrees of motion but also for the dynamical change of deformation, 
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because one might suppose that even if the trinuclear system is almost cold the 
small excitation energy present in the reaction will induce a /3-polarization of the 
fragments 0- 
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TABLES 


TABLE I. The tip distance d, the initial kinetic energy, E®, the final kinetic energy 
El and the asymptotic angle of the a, with point-like and with finite size Coulomb 
forces 


Splitting 

d(fm) 

E° a (MeV) 

El (MeV) 

Qf 

El (MeV) 

of 

w a 




Point-like forces 

Finite-size 

forces 


6 

2.71 

20.10 

80.83 

21.36 

77.86 

92 Kr+ 156 Nd 

7 

1.72 

15.84 

78.77 

16.91 

75.70 


8 

0.85 

11.25 

75.42 

12.10 

73.14 


6 

2.68 

19.87 

85.03 

20.96 

83.29 

104 Mo+ 144 Xe 

7 

1.70 

15.57 

83.77 

16.44 

82.08 


8 

0.82 

10.82 

81.63 

11.47 

80.32 


6 

2.85 

20.29 

88.17 

20.86 

86.99 

116 Pd+ 132 Sn 

7 

1.84 

15.86 

87.69 

16.31 

86.27 


8 

0.96 

11.53 

85.68 

11.19 

86.93 
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